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ABSTRACT 

Long term management of water resources in a catchment requires input from hydrological studies in the form of 
estimation or forecasting of the magnitude of hydrological variables. Such forecasts are useful in many ways, like warning 
of extreme floods or droughts, and help optimize the operation of systems, like reservoirs and power plants. This paper 
attempts to improve the original Soil Service Conservation Service Curve Number (SCS-CN) concept to perform 
long-term hydrologic analysis for a humid watershed, namely Hemavati, a tributary of River Cauvery in Karnataka State. 
Besides, the SCS-CN approach widely used as a simple method for predicting direct runoff volume for a given rainfall 
eventhas been evaluated for its applicability to humid catchment. The model inputs include daily rainfall, evaporation, 
and the model describes antecedent moisture (AM) by incorporating effect of antecedent rainfall without considering 
different antecedent moisture conditions (AMC). The model performs satisfactorily with efficiency of 87.51% and 88.71% 
in calibration and validation, respectively. The average relative error is of the order of 10%, which is tolerable. 
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1. INTRODUCTION 

Since the question of predicting and estimating the runoff resulting from precipitation and a quantitative 
understanding of the different processes involved in runoff generation are considered as fundamental and essential issues in 
hydrology, obtaining its quantitative and qualitative amount with a systemic approach is very important as it helps in 
decision making for the planning of water resources development and its management. Furthermore, runoff quantification 
gives indications on the opportunities to harvest rain water. 

In this study, original SCS-CN method has been modified by incorporating the effect of rainfall and used to 
compute direct surface runoff. The SCS-CN concept was selected because it was already adopted for various regions, land 
uses and climate conditions [1]. It provides an empirical relationship for estimating initial abstraction and runoff as a 
function of soil type and land use. Rainfall-runoff relationship can be visualized by the factors such as initial abstraction 
(la), direct runoff (Q), and actual retention (F). It considers the physiographic heterogeneity of the catchment (for example, 
topography, soil, and land use) to simulate the rainfall -runoff relationship at catchment level [2, 3]. The model has been 
widely used with success, providing consistently useful results [3, 4]. It was also evolved well beyond its original scope 
and it became an integral part of continuous simulation models [3, 5]. However, in spite of its widespread use, there is not 
an agreed methodology to estimate the CN values from measured rainfall-runoff data. Such a method would be important 
for two main purposes: (a) it would allow determination of CN values from measured rainfall -runoff data of local or 
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nearby similar watersheds when suitable data were available and (b) it would facilitate studies aiming at the extension of 
the SCS-CN method documentation for different, soil, land use, and climate conditions. CN varies as a function of the soil 
infiltration capacity and the land cover of the watershed, which are two essentially time invariant factors. The main 
difficulty is that the CN values calculated from measured rainfall-runoff data actually vary significantly from storm to 
storm on any watershed. Ponce and Hawkins (1996) reported that possible sources of this CN variability is the effect of the 
temporal and spatial variability of storm and watershed properties, the quality of the measured data, and the effect of 
antecedent rainfall and associated soil moisture [6]. Various sources of temporal variability, such as the effect of 
spatio-temporal rainfall intensity variability, the effect of antecedent rainfall, etc. make CN to be considered as a random 
variable. 

The methodology developed by the United States Department of Agriculture (USDA) has been documented in 
National Engineering Handbook, Section 4 (NEH-4) which has been revised several times (1956, 1964, 1965, 1969, 1972, 
1985, 1993, 2004) since its first publication in 1954. The SCS-CN method [2] was originally developed as a lumped model 
and up to this date it is still primarily used as a lumped model. In natural watersheds, however, spatial variability with 
regard to the soil-cover complex is inevitable (such spatial heterogeneity in the watershed could be considered temporally 
invariant). In this paper, an attempt has been made to develop an improved SCS-CN model to predict the catchment 
physics in a better way. 

Thus, the objectives of the present research are (i) a modification to the existing SCS-CN method to take care of 
the variability of daily CN values with an introduction of antecedent moisture factor due to the effect of antecedent rainfall; 
(ii) proposing a model to compute runoff components as well as other hydrologic components involved in the runoff 
production with a remark of dominancy/dormancy of each process. 

The proposed lumped conceptual rainfall-runoff model, designated as LCRR3 obviates the limitations of SCS-CN 
concept and is capable of simulating, other than direct surface runoff, the total stream flow and its components such as 
surface runoff, through flow, and base flow which is conceptualized to have two different moisture stores, i.e. soil moisture 
store and ground water store. This continuous simulation model considers a daily time step interval for analysis. A more 
general relation between I a -S including the effect of rainfall has also been introduced. 

2. EXISTING SCS-CN MODEL 

Soil Conversation Service (SCS) in 1975[2], based on multiple and numerous observations in famous fields and in 
different lands, has offered methods for estimating the runoff resultant from precipitation. The quantity of runoff (Q) is 
dependent on the precipitation (P) and the real holding (F). F is the difference between precipitation and runoff quantities. 
Moreover, a certain quantity of precipitation, at the start of rain storm, does not participate in runoff flow and is set aside 
for surface absorption, potholes, and penetrability capacity before the start of runoff, which is known as initial absorption 
(I a ).The existing SCS-CN method [4] is the sum of the following three equations [4]: 



P = I a+ F + Q 



(1) 



Q = F 

(p-i.) s 



(2) 
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h=M (3) 

where, P = total precipitation; I a = initial abstraction; F = cumulative infiltration; Q = direct runoff; S = potential 
maximum retention or infiltration; X = initial abstraction. The parameter S can be interpreted as post I a storage or the 
maximum depth of rainfall, excluding initial abstraction that can be potentially be extracted at a site[4].In order to simplify 
the equation and eliminate one variable, I a is fixed as0.2S. 

Combination of (1) - (3) leads to the following popular form of the SCS-CN method: 

Q= ( f - 7 «) 2 ,6= ^°^ (4,5) 
(P-I a +S) (P + 0.85) 

Here, P >I a and Q = otherwise. 

The above (5) shows that S is the only parameter that determines the volume of direct runoff. Quantity S is 
connected through a relation with an aspectless agent called CN. The quantity of CN is variable between and 100 and 
when CN becomes zero no runoff occurs, and when CN is equal to 100 all precipitation converts to surface runoff [7]. 
The retention parameter S is related to the value of CN by 

5 = 25.4f— -lo] (6) 
I CN J 

where S is in mm and CN is non-dimensional. Equation (6) shows that the CN decreases as S increases. 

3. FORMULATION OF SCS-CN-BASED LCRR3 MODEL 

This suggested lumped conceptual rainfall -runoff model, LCRR3, which is based on the modified concept of 
SCS-CN, attempts to simulate daily runoff from given daily rainfall and evaporation. The model is conceptualized to have 
two different moisture stores: soil moisture store and ground water store. The model consists of three major runoff 
components: (i) surface runoff, (ii) through flow and (iii) base flow. Out of these three, the surface runoff is computed 
based on the modified SCS-CN concept by incorporating antecedent moisture factor without considering different 
antecedent moisture conditions. The modified SCS-CN-based lumped model considers various hydrologic components 
involved in runoff generation mechanisms and takes into account the temporal variations of CN. The model formulation 
and various components of hydrologic cycles are described below. 

3.1 Initial Abstraction 

The initial abstraction I a is taken as a fraction of the possible retention in the soil and is computed as: 

I a ^ =X,S t , if t < 5 days (7) 

Here, X is the parameter to be optimised. I a(t) and S t are the daily initial abstraction and daily potential maximum 
water retention. Otherwise 



r a(t) - 



s, 



P.+S. 



, if t > 5 days (8) 



www.tjprc.org 



editor@tjprc.org 



22 



K. Geetha, S. K. Mishra, T. I. Eldho & Sangeeta Verma 



Here X and a are the coefficient and exponent of the I a which are to be optimised and P t is daily rainfall. 



3.2 Antecedent Rainfall 



The 5-day rainfall prior to the storm is considered as antecedent rainfall (ANTRF) and is computed as: 



ANTRF, = P (M) + P (f _ 2) + P (; _ 3) + P M + P M 



(9a) 



where t is the current day and P t is the rainfall of the respective day. 
Antecedent Moisture 

This model also considers the current space available for retention S t *for the first 5 -days assuming CN t = CN or 
S t *= S . The number of days exceeding 5, antecedent rainfall at any time't' (ANTRF t ) is computed from (9a). Using the 
antecedent rainfall value, the antecedent moisture (AM) amount can be computed as follows: 



Here, P is the parameter which is to be optimised. With the value of AM t , the current day possible water retention 
S t is modified as follows: 



Here S t * is the same as S t in (7), but corresponds to CN . The S, in (9c) is again modified by the 
evapotranspiration loss and drainage from the soil moisture zone and the daily input due to infiltration. 

3.3 Surface Runoff 

The amount of rainfall reaching the ground (P e ) after I a is responsible for initiating different kinds of paths of 
surface and subsurface flows. Replacing Q by RO t (daily surface runoff) for clarity in text, (4) can be re-written for daily 
runoff with time t as subscript as [4, 8-1 1]: 




(9b) 




(9c) 



RO t = 




P t -Ia(t)+S t 



or 



RO t = 




Pe(t)+S t 



(10) 



where 



P e(t) =Vt- l a(t) 



(11) 



I a(t ) =xs t 



(12) 
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S t = 2540 ^-254 (13) 
CN t 

whereP t =daily rainfall; I a(t) =daily interception; P e(t) =daily effective rainfall. The value of coefficient of initial 
abstraction X is assumed as 0.2 and also P e (,)^ else RO t = 0. Using the daily effective rainfall (P et ), the daily rainfall 
excess RO t can be computed by using (10) for the first 5-days of simulation, only if P exceeds I a , it is zero otherwise. 

3.4 Routing of Rainfall Excess 

When the number of days exceeds 5, to transform the surface runoff that is produced at the outlet of the basin, 
RO t (10) is routed using a single linear reservoir concept, as in (14). Then RO t is routed to the outlet of the basin using the 
single linear reservoir as below [4, 9-12]. 

SRO, =D x RO, + D, x RO, , + D 2 x SRO, , (14) 
where 

D =~^ d4a) 
2 + (\IK) 

D l = D (14b) 

D 2 ^-4^\ (14c) 
2 + {\IK) 

Here, SRO, is the routed surface runoff at the outlet of the catchment and K is the storage coefficient. 

3.5 Infiltration 

This amount of water reaching the ground after I a and not produced as surface runoff is assumed to infiltrate into 
the upper soil. F,_! is the previous day infiltration (mm) computed using water balance equation [4, 10- 11]: 

^(f-i) = *V-i) _ 7 «(r-i) _ R0 (t-1) (I 5 ) 

Here, ifP e(t) > 0, F,> 0. 

3.6 Evapotranspiration 

The amount of water lost as evapotranspiration can be computed by summing up the daily evaporation from the 
water bodies and transpiration from the soil zone in the watershed. 

3.6.1 Evaporation 

The daily evaporation EV, is computed as follows [4, 10-11]: 

EV t - PANC x EVP t (16) 

Where EVP, is the potential evaporation based on the field data and PANC is the Penmann coefficient, assumed 
as 0.8 for June-September and 0.6 for October-November, and 0.7 for February-May [10,11] 
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3.6.2 Transpiration 

Transpiration from the soil zone is considered as a function of water content available in the soil store above the 
wilting point of the soil [11]. It is computed as: 

TR t =K^{S abs -S t -e w ) (17) 

where Ki = coefficient of transpiration from soil zone, 9 W = wilting point of the soil, S abs = absolute maximum 
potential water retention, and S t = maximum possible water retention on t* day. The total actual evapotranspiration is taken 
as the sum of evaporation and transpiration as follows: 

ET t =EV t +TR t (18) 

3.7 Drainage 

The term drainage is used as the outflow from a linear reservoir only when the moisture content in the soil zone 
increases and exceeds the field capacity 9 f [11-13] as: 

DR t =K 2 x(S ahs -S t -0 f ) (19) 

where K 2 =subsoil drainage coefficient, 9 f = field capacity of the soil, and DR t is the drainage at any time 't'. 

3.8 Throughflow or Interflow 

The outflow from the unsaturated soil store is partitioned into two components: (i) subsurface flow in lateral 
direction and (ii) vertical percolation into ground water zone. The former component representing the through flow is taken 
as a fraction of the above drainage rate [11-13]: 

THR t =K 3 x DR t (20) 

Where THRt = through flow at time 't' and K 3 = unsaturated soil zone runoff coefficient. 

3.9 Percolation 

The outflow in the vertical direction from the unsaturated zone meets the ground water store due to the 
permeability of the soil. This percolated amount of water (PR t ) is considered as a part of drainage, and it is estimated as 
[11-13]: 

PR t ={l-K 3 )DR t (21) 

3.10 Deep Seepage 

The saturated store is considered as a non-linear reservoir and from this saturated store, outflow occurs at an 
exponential rate in the form of deep seepage DS t . This is modeled as follows [11-13]: 



DS,={w, -Vff 



(22) 



where DS t = deep seepage at any time 't'; v|/ t is the ground water content at any time 't'; v|/ f is the field capacity of 
the ground water store; and E = exponent of ground water store. 
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3.11 Base Flow and Deep Percolation 

The base flow of a watershed is the ground water release from a catchment in a stream. This active ground water 
flow which is also known as delayed flow can be modeled as outflow from a non-linear storage in the form of base flow 
(BF t ) [11-13]. 

The remaining amount of deep seepage which goes into aquifers in vertical direction is considered as a loss from 
saturated store and is taken as deep percolation DPr, [11-13]: 

BF t = BCOEF x DS t (23) 
D Pr, = (l - BCOEF )xDS t (24) 

Where BCOEF= ground water zone runoff coefficient, DPr t = deep percolation at any 't' and BCOEF = ground 
water zone runoff coefficient. 

3.12 Total Stream Flow 

The total stream flow (TRO t ) on a day 't' is obtained as the sum of the above three components, surface runoff, 
throughflow, and base flow [11, 12]. if t < 5days, 

TRO t = RO t + THR t + BF t (25a) 

and if t> 5 days, 

TRO, = SRO, + THR t + BF t (25b) 

3.13 Water Retention Budgeting 

The water balance in the soil and ground water store is worked out as follows [11, 12] 

de 



= F t -ET t -DR, 
dt 



(26a) 



^- = PR t - BF t - D Pr, (26b) 
dt 

dO dy/ 

Where and represent changes in water contents in soil moisture store and ground water store, 

dt dt 

respectively. This proposed model (designated as LCRR3) consists of 15 parameters: CN ,X, a, P, K, S a b S , 6f, 8w, Ki, K 2 , 
K3, v|/ fi v|/ t( i), BCOEF, and E. 

4. MODEL EFFICIENCY 

The Nash-Sutcliffe efficiency (NSE) was used to assess the SCS-CN model performance. NSE is a normalized 
statistic that determines the relative magnitude of the residual variance compared to the measured data variance [14]. 
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NSE indicates how well the plot of observed versus estimated data fits the 1:1 line. The efficiency of the model is 
computed using[14]: 



Efficiency = 



where 



x 100 (27) 



RV = f J {Q i -Q) (28) 
1=1 



2 



/v=£fe-a) (29) 

i=l 

Here, RV is the remaining variance; IV is the initial variance; Qj is the observed runoff for i th day; <2 ; is the 

computed runoff for i th day; n is the total number of observations; and Q i is the overall mean daily runoff. Efficiency is 
used for evaluating the model performance. Efficiency varies at the scale of to 100. 

4.1 Error Criteria 

The relative error (RE) is also computed to evaluate the deviation of simulated runoff from observed runoff 
[4, 10-12]: 

Relative error RE(%) = ^" x 100 (30) 

Qo 

Here, Q Q is the observed runoff and Q c is the computed (simulated) runoff. The higher RE value is indicative of 
greater deviation from the observed and vice versa. 

5. STUDY AREA AND DATA USED 

The study area selected is the Hemavati catchment, a tributary of Cauvery, originating in the Western Ghats in 
Mudgiri taluk of Chikmagalur district in Karnataka State, India (Figure 1). The drainage area of the catchment is 
600sq.km.For Hemavati catchment, daily data of three years i.e 1974-1976 are used for calibration and the remaining 
2 years (1977-1978) are used for validation of the model. Hydrologic data collected for this study consists of daily rainfall, 
evaporation and runoff data. 




Figure 1: Drainage Map of Catchment Hemavati 
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6. MODEL APPLICATION 

The modified SCS-CN model introduces the antecedent moisture factor with the effect of antecedent rainfall to 
estimate the daily water retention store and then the water retention store gets updated by taking into account of daily 
evapotranspiration, drainage from soil moisture store and infiltration to the soil moisture store. The optimal estimates of 
model parameters were obtained by using non-linear Marquardt algorithm coupled with trial-and-error [4]. 
The ranges/initial estimates are chosen appropriately. Table 1 shows the ranges and initial estimate of each parameter and 
also the optimised values of the parameters involved in the model formulation. The Nash &Sutcliffe efficiency along with 
the runoff coefficient is also presented in Table 1. It is seen that the catchment Hemavati shows runoff coefficient as 0.782 
and hence can be classified as a high runoff producing catchment. 

The model yields an efficiency of 87.51% in calibration, and 88.71%in validation, respectively. The higher 
efficiency reveals that the modified SCS-CN model is efficacious to high runoff producing Hemavati catchment. 

Table 2 presents annual values of rainfall, observed and simulated runoff and also computes error in percentages 
of runoff. It is seen that catchment Hemavati receives annual average rainfall of 2854.19mm, falling in humid region. 
While comparing the observed and simulated runoff (Table 2), it is seen that the model overestimates the runoff in the year 
1978-1979. The annual average relative error (RE) is also computed and it ranges from 0.21 to 25.78% with an average 
error of 10.41% (Figure 2). These values generally exhibit a satisfactory model performance. The RE value of -4.10% for 
1978-79 implies that the model overestimates runoff in this year. 

Table 1: Estimates of Parameters and EFFICIENCY of the Model LCRR3 



SI. No. 


Parameter 


Hemavati Catchment 


Range 


Initial 
Values 


Optimised 
Values 


1 


CN 


0-100 




45.00 


2 


X 


0-1.0 




0.001 


3 


a 


0-7.00 




5.00 


4 


P 


0.5-4.0 


1.0 


3.193 


5 


K 


0.1-2.0 


0.5 


0.93 


6 


K, 


0-1.0 


0.2 


0.03 


7 


K 2 


0-1.0 


0.01 


0.044 


8 


K 3 


0-1.0 




0.090 


9 


BCOEF 


0-1.0 




0.93 


10 


E 


0-1.0 


0.2 


0.487 


11 


Sabs 


200-900 


500 


835.78 


12 


W 


40-100 




60.00 


13 


e f 


200-500 


300 


319.41 


14 


% 


200-600 


300 


347.83 


15 


^t(l) 


40-300 




140.00 


16 


Efficiency 
(Calibration) 






87.51% 


17 


Efficiency 
(validation) 






88.71% 


18 


Runoff 
factor 






0.782 
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Table 2: Annual Rainfall, Observed Runoff, Simulated Runoff and 
Relative Error for Model LCRR3 Using Annual Data 



SI. 

No. 


Year 


Rainfall 

(mm) 


Observed 
Runoff 
(mm) 


Simulated 
Runoff 
(mm) 


Relative 
Error (%) 


1 


1974-75 


2937.50 


2552.53 


2128.36 


19.93 


2 


1975-76 


2650.89 


1717.96 


1682.78 


2.05 


3 


1976-77 


2676.33 


1894.46 


1890.46 


0.21 


4 


1977-78 


2941.87 


2936.71 


2179.53 


25.78 


5 


1978-79 


3064.35 


2061.74 


2146.24 


-4.10 


Average 


2854.19 


2232.68 


2005.47 


10.41 



Figure 2 shows the daily variations of estimated and observed runoff with respect to daily average rainfall for 
Hemavati catchment. While comparing, it is apparent that there is a good match between observed and simulated runoff. 

Model performs satisfactorily in the catchment except few peaks, where the computed runoff is lower than 
observed runoff. This is due to mainly the optimisation of the parameters by minimising the error based on a large number 
of other data points than peak. 



Hemavati (Jul* 1974-May1979| 
R£. 14 41H 



- Huioff |E( 
Rmfl (0> 

-Rainfall 



k k, ftl 




1 1S1 361 541 721 301 1081 1261 1441 1621 1M1 

Time in Day& 



Figure 2: Daily Variations of Rainfall, Observed Runoff (O), Estimated Runoff (E) with Average Relative 
Error (R.E) of Catchment Hemavati for Model LCRR3 Using Annual Data 




■ Rainfall (mm] 

■ Observed Runoff (mm) 
Simulate d Runoff (mm) 



Figure 2: Annual Estimates of Rainfall, Observed Runoff and Simulated Runoff 

Figure 3 presents the annual estimates of rainfall, observed runoff and simulated runoff in Hemavati catchment 
using 5-year data. Catchment receives maximum rainfall of 3064.35mm in the year 1978-79, and model computes the 
simulated runoff slightly more than the observed one in the same year. Except 1978-79, all other years, simulated runoff is 
found less than the observed. 

Table 3 presents the percent estimates of all hydrological processes involved in runoff production. This revised 
SCS-CN model computes the different components of runoff such as surface runoff, through flow and base flow as well as 
determines the dominancy/dormancy of the various processes. It is apparent from Table 3 that the dormant processes 
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(marked as *) are initial abstraction, deep percolation into aquifers and through flow. High amount of runoff is generated in 
Hemavati in the form of surface runoff and base flow. It is also observed that nearly 78% of rainfall is turned as runoff 
(observed) and the model computes as nearly 70% is transformed as runoff (simulated). This also verifies the model 
efficiency as the percent estimates of both observed and simulated runoff are matching. 

The advantages of the proposed method are more evident in Figure 2 and 3, where the good match between the 
observed and estimated runoff is demonstrated. As it can be clearly seen, satisfactory runoff predictions can be obtained by 
the proposed modified SCS-CN methodology. 

Table 3: Percent Estimates of Hydrological Components of Catchment Hemavati for Model LCRR3 



SI. No. 


Components 


Catchment 
Hemavati 


1. 


Rainfall (P t ) 


100 


2. 


Initial abstraction (I a(t) ) 


* 

1.54 


3. 


Effective rainfall (P e ( t )) 


98.46 


4. 


Infiltration (F t ) 


69.09 


5. 


Drainage ((DR t ) 


46.56 


6. 


Percolation (PR t ) 


42.36 


7. 


Deep seepage (DS t ) 


39.62 


8. 


Deep percolation (DPr t ) 


* 

2.93 


9. 


Surface runoff (SRO t ) 


29.37 


10. 


Throughflow (THR t ) 


* 

4.20 


11. 


Base flow (BF t ) 


36.69 


12. 


Simulated runoff 


70.26 


(TRO t ) 


13. 


Observed runoff (Q Q b S ) 


78.22 



Dormant Process 

This study estimates the annual water yield of various processes considered usually helpful in planning for 
utilization of resources and identification of dominant/dormant processes. Figures 4 and 5 present the percent estimates of 
runoff components as well as hydrologic processes involved in runoff generation mechanism like initial abstraction, 
infiltration, drainage, percolation, deep seepage, and deep percolation. 

It is apparent from Figure 4 that through flow is insignificant as compared to other runoff components such as 
direct surface runoff and base flow. Base flow is the major contributor to the total runoff. Figure 5 indicates that major 
amount of water gets infiltrated into the soil and is computed as nearly 69%. This study also reveals that the losses like 
initial abstraction and deep percolation into the aquifers are insignificant in the catchment. 



Percent estimates of Runoff 
components 

40 n 




a! SRO THR BF 



FLu n □ f f co n* p o n e n t :. 



Figure 3: Percent Estimates of Runoff Components of Catchment Hemavati 
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Figure 4: Percent Estimates of Hydrologic Processes of Catchment Hemavati 
7. CONCLUSIONS 

A lumped conceptual rainfall -runoff model was suggested with a modification to the existing SCS-CN original 
and applicability of this model was investigated in a humid catchment, River Hemavati, in Karnataka State, India to predict 
the behaviour of catchment. The model performance was tested with Nash-Sutcliffe efficiency. The following conclusions 
are drawn from the study: 

• The suggested long term hydrologic model modifies the value of daily potential water retention (S) by 
incorporating the effect of rainfall on antecedent moisture and thus a better relation between I a and S is developed. 

• The modified SCS-CN model obviates the limitations of the original SCS-CN model and is capable of simulating 
various runoff components other than direct surface runoff. Also it is capable of simulating various other 
processes involved in the runoff production. Dormancy/dominancy of each process in hydrologic cycle is also 
identified. 

• Observing the performance of the model application in catchment Hemavati, using annual data set, 
the investigation shows that the model with large number of parameters simulates the catchment response 
successfully. 

• The computed efficiencies in calibration and validation as well as average annual relative error indicate model is 
efficacious to high runoff producing catchment. 
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